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FOREWORD 


At  the  present  time,  thousands  and  thousands  of  people  around  the  worid  deal  with  ice,  snow  and 
pennafirost.  They  are  scientists,  educates,  engineos,  navigators,  meteorologists  and  odiers.  While  a 
small  fraction  of  these  people  contribute  to  die  knowledge  base  in  ice  physics,  all  of  them  use 
knowledge  from  it  firequently.  Moreover,  successful  aj^lied  research  is  based  upon  fundamental 
science— one  more  reason  for  ice  specialists  to  have  a  texdxxdc  on  ice  physics  on  their  desks. 

The  first  modem  ice  physics  text  was  Fletcher’s  book  on  The  Chemical  Physics  of  Ice  (1970). 
Fletcher’ s  book  is  in  typical  textbook  format:  it  is  reasonably  brief  and  easy  to  understand.  He  touched 
on  a  few  of  the  most  inqxntant  topics,  but  not  all  of  them. 

The  most  comprehensive  book  on  ice  physics  to  date  was  published  by  Hobbs  in  1974.  Hobbs 
considered  almost  all  of  the  basic  aspects  of  ice  as  understood  at  that  time.  Moreover,  he  described  and 
con^iared  several  (sometimes  opposing)  viewpoints.  This  fundamental  and  rather  large  (837  pages) 
book  is  commonly  known  as  the  “Ice  Bible”  by  specialists  in  the  field.  In  1974  and  1975,  two 
monographs  on  ice  were  produced  by  John  Glen.  These  were  briefly  and  clearly  written  and  reviewed 
almost  all  ice-related  subjects.  This  wwk  was  (and  in  some  respects  still  is)  a  magnificent  introduction 
to  ice. 

Finally,  in  1 98 1  Maeno  wrote  a  simple  popular  book  for  the  express  purpose  of  attracting  people’ s 
attention  to  the  subject. 

During  the  past  20  years,  a  signifreant  amount  of  new  experimental  and  theoretical  wwk  has 
appeared,  dramatically  changing  our  views  on  ice  physics.  As  a  result,  we  are  now  able  to  fcnmulate 
physical  laws  using  more  simple  and  direct  methods.  We  have  found  some  of  the  physical  models  used 
in  the  past  to  be  completely  wrong.  The  physics  of  ice  is  a  much  better  develcqied  subject  than  it  was 
20  years  ago. 

For  the  above  reasons,  we  feel  the  time  is  ripe  for  acontemporary  book  on  ice  physics,  incorporating 
the  known  and  proven  with  almost  20  years  wrath  of  material  not  covered  by  previous  works. 

I  have  tried  to  prepare  a  “readable”  book,  and  not  one  that  requires  the  reader  to  be  a  uniquely 
educated  person.  It  is  my  intent  to  present  the  material  in  such  a  way  that  any  reader  attracted  by  the 
title  Ice  Physics  will  be  able  to  comprehend  it.  This  is  quite  difficult  for  a  book  dedicated,  not  to  a 
particular  field  of  knowledge,  but  to  a  specific  material.  Indeed,  fra  ice  it  means  we  have  to  consider 
a  wide  variety  of  subjects,  including  quantum  chemistry,  solid  state  physics,  the  theory  of  elasticity, 
ionic  conductivity,  synchrotron  X-ray  topography,  crystal  growth,  the  physics  of  surfaces  and  more. 

The  primary  goal  is  to  produce  as  simple  a  book  as  possible  without  sacrificing  scientific  accuracy. 
Experimental  facts,  physical  ideas  and  theories  will  be  strongly  oiganized  and  bound  together 
cohesively.  The  reader  will  be  introduced  to  a  wide  variety  of  material  on  a  step  by  step  basis.  Then 
the  picture  will  be  whole. 

To  accelerate  materials  publication,  this  book  will  appear  first  in  the  form  of  a  series  of  joint 
CRREL-Dartmouth  reports,  later  to  be  published  in  CRREL’s  Monognqph  series,  on: 

1 .  The  structure  of  ordinary  ice 
Part  I:  “Ideal”  structure  of  ice.  Ice  crystal  lattice 
Part  n:  Real  structure  of  ice.  Defects 


2.  Electrical  properties  of  ice 

Parti:  Conductivity  and  dielectric  permittivity  of  ice 
Part  I!:  Advanced  tt^ics  and  new  physical  phenomena 

3.  Optical  properties 

4.  Electro-optical  effects  in  ice 

5.  Thermal  properties 

6.  Mechanical  properties  of  ice.  Elasticity  and  anelastic  relaxation.  Plastic  pr(^>efties.  Fracture  of  ice 

7.  Electromechanical  effects  in  ice 

8.  Surface  of  ice 

9.  Odier  forms  of  ice  and  their  properties 

10.  Ice  in  space 

1 1 .  Ice  research  laboratories 

The  reports  will  be  prepared  in  a  sequence  convenknt  to  the  author.  The  present  rqxxt  is  the  second  in  the 
series. 
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Structure  d  Ordinary  Ice  ^ 

Part  I :  Ideal  StriMture  of  Ice 

VICTOR  F.  PETRENKO 

INTRODUCTION 

Unlike  most  of  the  piccediiig  books  and  inonogn^)hs  on  ice  physics  (Fletcher  1970,  Gten  1974,  Hobbs  1974). 
this  report  will  be  concerned  only  with  the  structure  (tf  ofdinaiy  ice  I|,.  The  stnictiiies  of  numerous  other  fonns  of 
ice,  which  exist  either  at  very  high  pressures  or.  very  low  tai4)eratures,  will  be  considered  in  a  later  rqxxt  I  have 
chosen  such  an  approach  because  the  overwhelming  majority  of  people  interested  in  ice  deal  with  ordinary  hex¬ 
agonal  ice  I|,. 

The  other  essoitial  difference  in  the  description  of  ice  structure  is  the  separate  treatment  of  “ideal”  or  defect¬ 
less  structure  and  ‘Veal”  structure.  Such  a  division  reflects,  fffst,  the  important  role  that  drfects  play  in  nuitmals 
physics.  Second,  it  is  in  the  study  of  defects  (dislocations,  pdiU  defects,  grain  boundaries)  that  significant  ptxrgress 
has  been  achieved  since  the  last  monograph  on  ice  physics  was  published. 

A  new  topic  for  books  about  ice  will  be  a  section  descritring  the  structure  ofelectron  energy  bands  in  ice — this 
is  currently  a  ftonder  of  dieoretical  and  experimental  investigations  in  ice  physics.  Lattice  elastic  vibrations — 
phonons — ^vtill  be  discussed  in  the  report  devoted  to  the  thermal  properties  of  ice. 

CHEMICAL  CONTENT  OF  WATER  AND  ICE--IMPURrnES  AND  ISOTOPES 

Let  us  consider  briefly  some  basic  properties  of  liquid  water  prior  to  analyzing  ice  structure.  First,  water  is 
directly  present  in  many  of  the  instances  that  we  ate  ^ng  to  discuss:  for  example,  when  ice  is  growing  from  a 
melt  or  when  ice  is  melting  and  water  is  {xesent  at  the  ice  surface  in  the  form  of  a  thin  film. 

Second,  many  physical  prc^rerties  of  ice  and  water  are  similar.  This  is  attributable  to  water  being  a  mdecular 
liquid  and  ice  being  a  molecular  crystal.  The  fact  ttutt  water  molecules  in  these  substances  preserve  their  individ¬ 
uality  determines,  for  instance,  some  of  the  optical  and  thermal  properties.  In  dmse  cases  of  stMiie  major  differ¬ 
ence  in  the  behavior  of  water  and  ice,  this  gives  a  key  to  understanding  the  mechanism  of  the  phenomenon,  since 
the  structural  distinctions  between  ice  and  water  can  be  formulated  quite  accurately. 

There  is  no  point  in  proving  to  the  reader  that  water  consists  of  oxygen  and  hydrogen,  since  we  acquire  this 
knowledge  in  high  school.  One  of  the  most  simple  and  accessible  ways  to  determine  the  chemical  content  of 
water  is  via  electrolysis.  I  believe  that  many  readers  successfully  performed  this  simple  eiqieriment  back  in  dwir 
childhood.  By  measuring  the  volumes  of  gases  generated  in  the  electrolysis  process 

2H2O  — >  2H2  +  O2  (1) 

and  using  Avogadro’s  law,  we  can  establish  the  chemical  formula  of  water:  H2O.  One  should  keep  in  mind, 
however,  that  the  path  to  understanding  this  at-frrst-glance  sinqile  result  was  tiresome  and  difficult  Major  omtri- 
butions  to  this  research  were  made  by  such  famous  scientists  as  Cavendish,  Lavoisier,  Dalton  and  Avogadro. 
Water  is  a  strong  and  fairiy  versatile  solvent  Therefore,  water  found  in  nature  usually  contains  substantial  amounts 
of  various  salts,  gases  and  organic  compounds  in  dissolved  form.  Owing  to  the  high  solubility  of  oxygen  in 
water,  thousands  of  diverse  animals  live  in  seas,  rivers  and  lakes. 


Talik  1.  Flijricalprapertktof  ordlnrjr  (B]0)  mmI  heavy  (DjO)  water  (after 
hfakakovlWS). 


Parameter 

HrO 

Ocasity 

iceW0°C 

0.917  t/oa? 

1.017  g/enP 

liquidwZO^ 

0.9982  gfem^ 

1.1053  gfcaP 

Vitcocity  at  20*C 

1.005  c 

1.251c 

MeMag  poial 

273.16  K(0%) 

276.97  K  (3.815%) 

BoU^  pout 

373.16  K  (100%) 

37439  K  (101.43%) 

Critical  tpnpoature 

647.3  K  (374.15%) 

643.9  K  (370.7%) 

Critical  pRMure 

22.06  MPa 

21.67  MPa 

Critical  deaaty 

0.322  gtaa^ 

0.356  g/ciiP 

Laieat  beat  of  ftwiaa 

332.4  l/g 

316.6  J/g 

Lateai  hew  of  vaporatioa 

2256.2  J/g 

2070.9  J/g 

latent  beat  of  aiWiiiiation  of  ice  at  0*C 

2833.91/8 

2639.6  J/g 

Specific  heat  capacity  (cj 

icettO*G 

2.038  l/gK 

2.202  J/g  K 

liquid  at  0°C 

4.186  i/gK 

4.23  J/g  K 

vapor  at  0°C 

1.905  i^K 

1.68  J/g  K 

Dielectric  penaiBivity 

iceW-10*C 

95 

92 

liquid  at  2S'*C 

78.54 

78.25 

Thennal  conductivity 

ice  at  0^ 

234.6  mW/mK 

— 

liquid  at  0*C 

560.9  mW/m  K 

559.5  mW/in  K 

liquid  at  100°C 

678.6  mW/m  K 

644.4  mW/m  K 

vapor  at  100% 

25.1  mW/rii  K 

— 

Adiabatic  compteatfinlity  W  20°C 

4JS5  X  lfr->7  mW 

4.70xl0-'lm5/N 

Dielectric  lelaxatioa  time 

ice  W -10% 

5  X  10-5  s 

9.1  X  10-5  s 

liquid  at  25% 

9.22x10-'^$ 

11.89xl0-'5s 

Molecular  magnetic  susceptibility  at  20% 

-12.972  xl0-« 

-12.948  X  10< 

Surface  tension  of  liquid  water 

atO% 

74.64  X  10-5  N/m 

— 

at  20% 

72.75  X  10-J  N/in 

72.57  X  10-5  N/m 

attOO"C 

58.89  X  10-5  N/m 

58.85  X  10-5  NAn 

Re&active  index  at  20°C 

1.333 

1.328 

The  chemical  fonnula  H2O  doesn’t  exhaust  the  con- 
tent  of  “absdutely”  chemically  pure  water  because  of  the  i8.oes  V 

presence  of  hydrogen  and  oxygen  isotopes.  *§  \  / 

Three  hydrogen  isotopes  are  known:  •H-oidinaiy  hy-  w  '  \  / 

drogen;^  or  D,  deuterium;  and  or  T.lritiunL  In  natu-  1  ie.oei  ■  \  / 

ral  wato’  thoe  are  IS  deutmum  atoms  per  10^  hydrogen  I  \  / 

atoms.  Tritium  is  ^-radioactive.  It  disintegrates,  electrons  3  \  / 

and  neutrinos  being  emitted  from  die  nucleus,  and  trans-  *  18.017  - 

forms  into  ^He  with  a  half-life  of  12.43  years.  Minute  . 

amounts  of  tritium  water,  (melting  temperature  4.S°C),  -4  -2  0  a  4  e  a  10  12 

^  TBJPERATURE  (C^\ 

can  be  either  synthesized  or  found  in  small  quantities  near 

some  radioactive  sources,  while  so-called  heavy  water  Figure  1.  Temperature  dependence  of  a  molar 

(riiairdyintheformofD20)isanintegralpart(althougha  ^ 

small  one,  0.015%)  of  the  world.  At  present  it  is  widely  (c0er  Malenkov  1988). 

used  in  various  industries  and  scientific  research. 

There  are  also  six  known  oxygen  isotopes:  '^o,  I'^O,  1*0  and  Among  them,  only  three  are 

stable:  and  **0.  For  every  10*  atoms  of  **0,  there  are  20  atoms  of  **0  and  4  atoms  of  •‘^O. 

Some  basic  physical  properties  of  ordinary  (H2O)  and  heavy  (D2O)  water  are  summarized  in  Table  1  (Malen¬ 
kov  1988). 

One  of  the  surprising  and  peculiar  properties  of  wato-  is  the  presence  of  a  density  maximum  at  4°C  (Ftg.  1). 
Because  of  this  maximum,  water  freezes  from  the  surface  rather  than  from  bottom.  It  is  hard  to  imagine  what 


Figure  1.  Temperature  dependence  of  a  molar 
volume  of  water  in  the  region  of  the  density 
maximum  {after  Malenkov  1988). 
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ecological  catastrophe  might  result  if  tbm  woe  no  such  “otinof'’ peculiarity  (rf  water.  Despite  the  widely  accept- 
ed  opinion  that  such  behavior  is  unique,  similar  peculiarities  in  the  density  of  a  liquid  phase  are  found  in  the  case 
of  In2Te3  (Stillingo’  1990). 

Analogous  density  maxima  are  obsoved  in  the  case  of  DjO  at  7=  1 1.2°C  and  T2O  at  7=  13.4°C  (StUlingo' 
1990). 

Anotho-  unusual  pt(^)erty  of  water — decrease  of  specific  volume  at  melting  (8.3% ) — is  also  vital  for  ecolog¬ 
ical  equilibrium.  Because  of  this  ice  floats  on  the  wato'  surface.  This  propoty  is  not  unicpie  either,  it  is  also 
observed  in  Ge,  Si  and  Bi. 

The  structure  of  ordinary  wam  will  be  considered  in  the  Position  of  Oxygen  Atoms  section,  along  with  the 
structure  of  ice,  and  a  comparison  analysis  will  be  performed. 


THE  WATER  MOLECULE 


One  of  my  colleagues,  an  expert  in  physical  chemistry,  said  once  that  he  had  spent  25  years  studying  the 
properties  of  the  water  molecule.  He  realized  that  his  whole  life  was  not  enough  to  answo-  all  the  questions  about 
this  molecule,  consisting  of  only  three  atoms.  Taking  into  account  his  expoience,  we  shall  abide  by  consider¬ 
ations  of  reasonable  sufficiency  for  the  test  of  the  material. 

We  shall  begin  our  acquaintance  with  the  water  molecule  with  a  considera¬ 
tion  of  its  external  charactoistics — its  shape  and  dimensions.  The  main  sources 
of  the  information  discussed  below  are  X-ray  diffraction,  c^cal  and  electron 
scattering  spectra  of  water  vapor,  measurements  of  the  heat  capacity  of  water 
vapor  and,  finally,  thetnetical  calculations. 

The  first  thing  that  strikes  and  excites  the  imagination  when  one  looks  at  a 
water  molecule  is  its  triangular  stuqie.  The  arrangement  of  oxygen  and  hydro¬ 
gen  nuclei  in  this  molecule  is  shown  in  Figure  2.  At  fust  glance,  such  a  strange 
shape  with  low  symmetry  (compared  to  a  linear  atomic  arrangement)  is  not  jus¬ 
tified;  it  violates  the  principles  of  simplicity  and  ratiwiaiity  that  we  admire  so 
often  in  the  world  around  us. 
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Figure  2.  Arrangement  of 
oxygen  and  hydrogen  nuclei 
in  water  molecule. 


Fortunately,  the  schematic  in  Figure  2  doesn’t  reflect  all  the  features  of  the  H2O  molecule’s  structure,  which 
are  determined  by  the  exterior  contours  of  electron  clouds  of  the  molecule.  As  we  shall  see  further,  the  sh^  of 
these  electron  clouds  is  symmetrical  to  a  very  large  extent  and  is  able  to  satisfy  our  unconscious  aspiration  to 


harmony. 

The  arrangement  of  atoms  in  the  vertices  of  a 
triangle,  represented  in  Figure  2,  is  not  the  only 
possibility.  Alternative  configurations  of  the  H2O 
molecule  appear  when  the  molecule  is  excited  (for 
instance,  by  ultraviolet  light)  and  when  the  mole¬ 
cule  undergoes  transition  from  the  ground  state 
(with  minimal  energy)  to  one  of  the  excited  states 
(Claydon  et  al.  1971). 

Table  2  contains  values  of  OH  bond  lengths 
and  the  angles  between  them,  calculated  for  sev¬ 
eral  energy  states  of  the  H2O  molecule.  At  this 
point  we  are  not  going  to  focus  on  the  meaning  of 
the  chemical  symbols  denoting  different  states.  For 
us  it  is  only  essential  to  know  that  the  symbol  rep¬ 
resenting  the  ground  state,  with  which  we  shall 
mainly  deal,  is  situated  in  the  upper  row  of  this 
table,  and  other  states  differ  from  the  ground  state 
by  energies  of  several  electronvolts.  As  can  be  seen 


Tabk  2.  EquiUbrium  geomctiy  of  the  lowest  dectronk  states 
of  H2O  with  respect  to  symmdrical  stretchfaig  and  bending 
coordinates  (after  Claydon  etaL  1971). 

VerticcU 


State 

Valence  MO 
configuration 

Bond 

length 

(A) 

Bond 

angle 

(deg.) 

Energy 

minimum 

(tV) 

excitation 

energy 

(eV) 

Singlels 

X'A, 

3a, 3  lb, 2 

1.03o 

104 

0.00 

— 

A<B, 

3a,2|b,4a, 

I.II5 

180 

5.79 

7.43 

'A2 

3a,2|b,22b2 

l.llo 

76 

6.93 

8.03 

B'A, 

3a,  lb, 24a, 

111, 

180 

5.79 

10.85 

•Bj 

3a,lb,22b2 

lIOo 

122 

10.27 

11.06 

Triplets 

3B, 

3a,2lb,4a, 

I.II7 

180 

4.76 

6.20 

3A2 

3a,2lb,2b2 

I.II7 

80 

6.55 

7.49 

3a,  lb, 24a, 

l.llo 

180 

4.76 

8.67 

3Bj 

3a,tb,22b2 

_Lil5_ 

112 

8.90 

9.43 
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from  Table  2,  among  excited  states  dioe  are  states 
with  a  linear  configuration  (angle  H-O-H  is  180°).  To 
understand  the  factors  that  govern  the  shape  of  the 
molecule,  one  has  to  consider  the  configuration  of 
electron  clouds.  Remember,  the  exact  solution  of 
Schrodinger’s  equation  for  such  a  con^licated,  many- 
body  system,  like  the  water  molecule,  doesn’t  exist. 
Several  highly  effective,  con^xitational,  semi-eiT^}ir- 
ical  techniques  have  been  developed  that  allow  us  to 
calculate  with  good  accuracy  the  reciprocal  arrange¬ 
ment  of  the  oxygen  and  hydrogen  nuclei  and  the  elec¬ 
tron  density  distribution,  as  well  as  energy  levels  in 
this  molecule. 

The  electron  cloud  of  the  hydrogen  atom  consists  of  a  single  election,  occupying  a  spherically  symmetrical  1  s 
orbital.  In  the  oxygen  atom,  two  inner  electrons  are  located  at  the  Is  orbital,  two  at  the  2s  orbital  (also  spherically 
symmetrical),  and  the  other  four  electrons  at  the  external  2p  orbitals.  Electrons  of  these  2p  orbitals  are  divided 
among  three  states  (2p2  has  two  electrons;  2py  and  2p^  have  one  election  each). 

According  to  modem  concepts,  when  two  hydrogen  atoms  with  their  electrons  are  attached  to  an  oxygen 
atom,  eight  electrons  take  part  in  the  formation  of  external  molecular  orbitals:  two  from  the  hydrogen  and  six 
from  the  2s  and  2p  orbitals  of  oxygen.  This  process  is  thus  called  s-p  hybridization.  The  water  molecule  electron 
cloud,  consisting  of  four  orbitals,  each  having  two  electrons,  results  from  this  hybridization  (Fig.  3). 

The  exterior  shape  of  the  molecule,  depicted  in  this  figure,  is  determined  by  the  electron  density  in  the  region 
around  the  oxygen  and  hydrogen  nuclei.  The  molecule  resembles  a  four-homed  snail,  the  horns  of  which  point 
approximately  in  the  directions  of  the  vertices  of  a  regular  tetrahedron.  The  H-O-H  angle  in  the  water  molecule  is 
equal  to  104.523°,  while  the  conesponding  angle  in  a  regular  tetrahedron  is  109.467°.  The  distance  OH  in 
equilibrium  is  0.958  A. 

The  two  molecular  orbitals  in  which  the  hydrogen  nuclei  are  located  are  called  “bonding  orbitals,”  while  the 
two  others  are  called  "lone-pair  orbitals.” 

The  electron  density  contour  map  for  valence  electrons  in  the  plane  H-O-H  of  a  water  molecule  is  shown  in 
Figure  4  (Bader  and  Jones  1963). 

The  nature  of  the  bond  between  oxygen  and  hydrogen  in  the  water  molecule  is  of  particular  interest.  In  the 
first  place,  since  they  share  electrons,  or  “nationalize”  them,  one  can  say  that  this  is  covalent  bonding.  However, 
since  a  part  of  the  electron  cloud  is  “deferred”  from  the  hydrogen  atom  to  the  oxygen  atom,  they  assume  opposite 
electric  charges  (+  at  the  hydrogen,  -  at  the  oxygen).  The  electrostatic  atbaction  arising  in  this  case  is  typical  for 
ionic  bonding.  Therefore,  one  can  talk  about  mixed  ionic-covalent  bonding  in  the  case  of  the  water  molecule. 
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Figure  4.  Electron  density  contours  in 
the  water  molecule  (after  Bader  and 
Jones  1963). 


From  our  knowledge  of  electric  charge  distribud(»  in  the  water  molecule,  the  dipcde  moment  can  be  calculat¬ 
ed  and  conqKued  to  the  value  that  can  be  derived  frcxn  data  on  the  dielectric  permittivity  of  water  vapor.  Fairiy 
good  conformity  exists  b^een  theoretical  and  expoimental  values 

Po  =  (6.1±0.02)xl0-30Cm.  (2) 

Owing  to  the  presence  of  more  than  two  point  charges  in  a  water  molecule,  it  possesses  multipole  moments  of 
higher  orders  than  dipole:  quadrupole  and  octupole.  Though  we  are  not  grang  to  employ  these  quantities,  the 
reader  can  find  detailed  information  about  them  in  the  paper  by  Glaeser  and  Coulson  ( 1%5),  for  instance.  Let  us 
alsonotethatsiiKethewatermoleculedoesnothaveanuncoinpensatedmagneticmomentinthegroundstate.it 
is  diamagnetic  (see  Table  1). 

Because  of  large  bonding  energies  of  hydrogen  atoms  and  electrons  in  the  water  molecule,  water  is  a  very  sta¬ 
ble  chemical  compound.  The  energy  of  formation  of  the  water  molecule  fiom  an  O  and  2H  abMns  at  0  K  is  9.S1 
eV.  The  energy  required  to  remove  a  hydrogen  atom  fiom  the  oxygen  atom  is  S.  1 1  eV  for  the  first  atrnn,  and  4.40 
eV  for  the  second  atom.  Photodissociation  of  watn-  molecules  (in  water  vapor)  begins  at  photon  energies  2  7.S1 
eV  (A,  S 165  nm). 

At  this  point,  we  finish  our  short  review  of  this  imptxtant  molecule.  Properties  of  the  H2O  molecule  associat¬ 
ed  with  rotations  and  oscillations  will  be  considered  in  the  rqxirt  on  the  optica]  ptiperties  of  ice. 

THE  HYDROGEN  BOND 

The  water  molecule,  described  in  the  previous  section,  is  the  block  that  we  shall  use  to  built  the  crystal  lattice 
of  ice.  As  can  be  seen  from  Figure  3,  the  'Mata  molecule  ideally  fits  the  role  of  construction  element  for  building 
an  ice  lattice.  This  molecule  already  has  four  “hands” — i.e.,  four  molecular  orbitals— stretched  towards  its  neigh¬ 
bors.  Two  of  those  hands  have  a  positive  charge,  owing  to  the  presence  of  protons,  and  the  other  two  have  a 
iwgative  charge.  This  allows  the  considoable  electrostatic  forces  acting  between  molecules  to  be  engaged  when 
molecules  associate  to  form  a  crystal.  Each  individual  molecule  forms  bonds  with  four  neighboring  molecules, 
as  is  illustrated  in  Figure  S. 

Two  adjacent  molecules  appear  to  be  connected  via  the  hydrogen  atom.  This  type  of  bond  is  usually  denoted 
O-H  -O  and  is  called  a  hydrogen  bond.  Other  atoms,  as  well  as  oxygen,  can  form  hydrogen  bonds.  Thus,  the 
general  symbol  for  the  hydrogen  bond  is:  A-H -  B.  Atexns  A  and  B  must  be  electronegative  and  there  must  be  a 
strong  covalent  bond  between  one  of  them  and  the  hydrogen  atom. 


region  of  high  / 

electron  probability- 


Figure  5.  Formation  of  four  hydrogen  bonds  between  a  water  molecule 
and  its  four  nearest  neighbors. 
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Figure  6.  H-bond  formation  by  H2O  and  HF  revealed  by  melting 
and  boiling  point  trends. 


To  the  veiy  first  approximation,  hydrogen  bonding  arises  from  the  electrostatic  interaction  between  an  elec¬ 
tronegative  atom  and  a  proton.  As  a  result,  the  distance  between  electronegative  atoms  becomes  approximately 
equal  to  the  sum  of  Van  der  Waals  atomic  radii;  in  other  words,  they  act  as  if  there  were  no  hydrogen  atom 
between  them  at  all.  For  example,  in  the  case  of  ice,  the  distance  between  two  adjacent  oxygen  atoms  =  2 .76 
A,  which  is  rather  close  to  the  sum  of  Van  der  Waals  radii  for  oxygen;  2  x  1 .4  A  =  2.8  A. 

The  energy  of  the  hydrogen  bond  is  1  to  1 .5  orders  of  magnitude  lowor  than  the  energy  of  a  typical  chemical 
(covalent)  bond,  but  at  the  same  time  it  is  2  to  3  orders  of  magnitude  higher  than  the  conesponding  energy  of 
non-valent  Van  der  Waals  interaction.  An  extremely  strong  hydrogen  bond  develops  between  HF  molecules, 
which  can  associate  to  form  polymer  structures  H2F2,  H3F3,..H^(.  Let  us  note  also  that  a  hydrogen  bond  can  be 
not  only  intermolecular,  but  intramolecular  as  well.  Such  intramolecular  bonds  are  one  of  the  major  factors  stabi¬ 
lizing  globular  proteins.  They  also  considerably  affect  the  properties  of  wood  and  paper  consisting  of  cellulose 
fibers,  and  are  responsible  for  the  unique  structure  of  nucleic  acids. 

Because  of  the  presence  of  a  branched  and  rather  dense  net  of  hydrogen  bonds,  ice  is  a  perfect  model  material 
for  studying  a  variety  of  processes  (proton  transfer,  for  instance)  that  occur  in  more  complex  organic  compounds. 

As  one  of  the  parameters  for  characterizing  the  strength  of  hydrogen  bonds  between  molecules,  the  values  of 
melting  temperature  and  boiling  temperature  T^,  can  be  taken.  In  Figure  6,  values  fra-  T„f  and  for  hydrates 
of  bivalent  and  monovalent  compounds  are  plotted.  As  is  clear  from  this  figure,  water  possesses  an  exceptional 
position  here,  having  the  highest  melting  and  boiling  temperatures.  In  common  ice  (I),),  the  energy  per  hydrogen 
bond  is  0.306  eV  (Whalley  1957,  1976).  This  quantity  is  calculated  as  half  the  energy  required  to  move  a  mole¬ 
cule  from  ice  into  the  vapor  (the  half  arises  horn  the  number  of  hydrogen  bonds  contributed  by  an  individual 
water  molecule  to  the  overall  structure). 


PHASE  DUGRAM  OF  WATER  AND  ICE 

As  mentioned  previously,  here  we  shall  deal  mainly  with  the  structure  of  ordinary  ice,  denoted  as  I|,.  The 
subscript  h  indicates  that  the  structure  has  a  hexagonal  symmetry.  Water  is  an  extraordinary  substance  in  many 
respects,  including  the  number  of  solid  phases — at  present,  14  of  them  are  known.  Therefore,  we  will  character¬ 
ize  the  temperature  and  pressure  ranges  in  which  ice  1|,  exists,  and  the  conditions  under  which  other  modifica¬ 
tions  of  ice  are  found. 

The  phase  diagram  of  ice  in  the  pressure  range  0  to  2.5  GPa  and  the  temperature  range  -2(X)  to  -100°C  is 
represented  in  Figure  7  (Whalley  1976).  We  shall  give  a  detailed  treatment  of  the  structure  and  properties  of  high 
pressure  phases  in  anotiier  report.  Here,  we  shall  only  point  out  that  all  ice  phases  shown  in  Figure  7  are  charac¬ 
terized  by  the  water  molecule  retaining  its  individuality  and  yet  being  linked  by  hydrogen  bonds  to  four  neigh¬ 
bors.  The  closest  (in  terms  of  the  structure)  “relative”  of  ice  I),  is  cubic  ice  4.  Cubic  ice  4  can  he  obtained  by 
water  vrqxir  sublimation  onto  a  substrate  in  the  temperature  interval  from  -80  to  -150°C  (Honjo  et  al.  1956).  Ice 
4  can  also  be  obtained  when  pressure  is  removed  from  some  of  the  high  pressure  phases.  Thus,  Bertie  et  al. 
( 1963)  discovered  that  ice  4  can  be  produced  in  large  quantities  “by  first  preparing  ices  II,  in  or  V  at  high  pres¬ 
sures,  cooling  them  to  liquid  nitrogen  temperatures,  reducing  the  pressure  to  atmospheric,  and  then  warming." 
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Another  ciystallogr^)hic  phase  of  ice,  named  ice  K,  forms 
during  the  slow  cooling  of  KOH-doped  ice  I),  below  72  K 
(Tajima  et  al.  1982).  Finally,  at  enomnus  pressures  P  ^  44 
GPa,  so-called  symmetrical  ice  (ice  X)  arises,  in  which  pro¬ 
tons  occupy  symmetrical  positions  with  respect  to  the  clos¬ 
est  oxygen  atoms.  The  relative  positions  of  ices  X  (symm^- 
rical-bond),  VUI,  Vn,  VI  and  the  liquid  phase  are  shown  on 
the  P-T  diagram  in  Hgure  8,  adopted  fiom  a  theoretical  pa¬ 
per  of  Stillinger  and  Schweizer  (1983).  In  fact,  the  division 
of  individual  molecules  vanishes  in  ice  X,  because,  owing 
to  the  symmetrical  position  of  the  hydrogen  with  respect  to 
the  oxygen,  it  is  impossible  to  tell  to  which  molecule  a  par¬ 
ticular  hydrogen  atom  belongs. 

Among  the  12  crystallographic  forms  of  ice  (ice  4,  I|,,  n, 
...,  XI),  only  5  (H,  VII-XI)  are  real  crystals,  in  which  a  trans¬ 
lational  symmetry  (the  pattern  is  repeated  when  one  goes 
from  one  unit  cell  to  another)  holds  for  both  oxygen  and 
hydrogen  atoms.  In  other  modiflcations  of  ice,  the  hydro¬ 
gen  subsystem  is  disordered.  We  shall  consider  this  “disor¬ 
der”  in  detail  in  the  Arrangement  of  Protons  and  Configu¬ 
rational  Entropy  of  Ice  section. 

The  other  two  solid  phases  (which  extend  the  sequence 
to  14)  are  ariKMphous  forms  of  ice.  The  first  one,  low-densi¬ 
ty  ariHxphous  ice,  can  be  produced  when  vapor  is  condensed 
onto  a  substrate  at  a  temperature  of  about  -1  SOX.  The  (ten¬ 
sity  of  this  phase  is  0.94  Mg  nr^.  The  second  kind  of  amor¬ 
phous  ice,  high-density  amcxphous  ice,  forms  at  T  =  77  K 
and  pressure  P  2  1  GPa  Its  density  (renormalized  to  a  pres¬ 
sure  of  1  atm)  is  1.17  Mg  m-3  (Whalley  et  al.  1987). 

We  also  will  be  interested  in  the  behavior  of  ice,  water 
and  water  vapor  in  the  vicinity  of  the  so-called  triple  point 
(+0.01X,  610  Pa),  at  which  all  three  phases  exist  in  equi¬ 
librium.  Figure  9,  adapted  from  Maeno  (1981),  shows  the 
phase  diagram  of  water  in  the  vicinity  of  this  triple  point 
and  phase  transformations. 


ARRANGEMENT  OF  OXYGEN  ATOMS 

I  have  no  doubt  that  readers  of  this  report  have  repeated¬ 
ly  noted  the  hexagonal  shape  of  most  snow  crystals  (Rg. 
10).  It  is  fairly  easy  to  deduce  that  this  shape  arises  from  the 
symmetry  of  ordinary  ice  crystals  at  the  molecular  level.  How 
does  such  hexagonal  symmetry  arise? 

To  understand  this,  let  us  consider  a  dense  packing  of 
identical,  ideal  balls  on  a  level  surface  (Fig.  1 1). 

After  the  first  layer  (let  us  denote  it  as  A)  and  the  second 
layer  (B)  are  packed,  there  are  two  alternatives  left  for  the 
third  layer.  First,  the  centers  of  the  balls  in  this  third  layer 
can  be  either  right  above  the  centers  of  the  balls  in  the  first 
layer,  i.e.,  in  positions  A,  or,  second,  in  positions  C. 


PRESSURE  (kBar) 


Figure  7.  Phase  diagram  of  ice  in  the  pressure 
range  0  to  2.5  GPa  (0to2S  kBar )  and  the  temper¬ 
ature  range  -200 to  -I00°C  (after  Whalley  1976). 


Figured.  Very  high-pressure  region 
of  the  H2O  phase  diagram  (after 
Stillinger  and  Schweizer  1983). 


PRESSURE  (Mbai) 

Figure  9.  Low-pressure  region  of  the  H2O 
phase  diagram  (after  Maeno  1981 ). 
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Figure  10.  Hexagonal  shape  of  snow  flakes. 


Figure  1 1.  Dense  packing  cf  identical,  ideal  balls  on 
a  level  surface. 
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If  consecutive  layers  of  balls  are  packed  in  the  order ...  ABABAB...,  then  the  structure  known  as  hexagonal 
close-packed  will  result.  If  we  now  replace  the  balls  with  dumbbells,  consisting  of  two  oxygen  atoms  connected 
by  a  hydrogen  bond  normal  to  the  plane  of  the  layers,  then  we  obtain  an  almost  exact  arrangentent  of  oxygen 
atoms  in  the  lattice  of  ice  Ii,.  Such  a  structure  in  two  projections  is  shown  in  Figure  12. 

The  arrangement  of  oxygen  atoms  in  the  ice  crystal  lattice  was  derived  from  numerous  data  from  X-ray  dif¬ 
fraction.  The  major  source  of  scattering  of  X-rays  in  a  crystal  is  electrons.  Therefore,  the  oxygen  atom,  surround¬ 
ed  by  10  electrons  (eight  of  its  own  and  two  from  hydrogen)  is  quite  an  efficient  source  of  scattering,  from  which 
fairly  good  X-ray  photogrtq>hs  may  be  obtained.  The  situation  is  quite  unsatisfactory  fw  the  hydrogen 
atom,  which  is  virtually  “striiq>ed”  of  its  electrons.  This  is  why  the  arrangement  of  hydrogen 
atoms  in  ice  was  inferred  mainly  from  otha  sources— neutron  diffraction  and  NMR. 

When  comparing  the  ice  Ih  structure  to  the  structure  of  close-packed  balls,  one 
should  take  into  account  that  the  ice  lattice  is  slightly  squashed  along  the  c- 


b.  View  is  along  the  hexad  axis  c. 


a.  View  in  a  direction  perpendicular  to  hexad  axis  c; 
minittutm  vectors  of  translation  are  indicated  as  c  arul  a 

Figure  12.  Positions  of  oxygen  atoms  in  the  ice  /*  crystal.  The  shown  size  of  atoms  is  not  in  real  proportion  to 
the  inter  atomic  distance.  Solid  lines  represent  hydrogen  bonds. 
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TiMr  1  Tf  niy  iHffrnrHtrn  mnmimnritfn  nf  thr  unit  rrl  pf  ■■iff n  nf  Iri  ■ 

ftmctkin  of  tempenrtuK. 


Ttmperatun  og  Cg  Tempenaitn  og  Cg 


(°C) 

(k) 

(k) 

(/) 

(k) 

From  LaPtaca  and  Post  (1960) 
-10  4.5190  7.3616 

1.6291 

Aoffl  Brill  and  Uppe  (1967) 

-80  4.5108  7.3438 

1.62804 

-20 

4.5169 

7.3570 

1.6288 

-100 

4.5062 

7.3360 

1.62798 

-40 

4.5128 

7.3500 

1.6287 

-120 

4.5025 

7.3298 

1.62794 

-60 

4.5088 

7.3440 

1.6289 

-140 

4.5000 

7.3261 

1.62802 

-80 

4.5052 

7.3388 

1.6290 

-160 

4.4987 

7.3240 

1.62803 

-100 

4.5021 

7.3344 

1.6291 

-180 

4.4979 

7.3225 

1.62798 

-120 

4.5001 

7.3304 

1.6289 

-200 

4.4974 

7.3215 

1.62794 

-140 

4.4990 

7.3267 

1.6285 

-220 

4.4971 

7.3208 

1.62789 

-160 

4.4961 

7.3233 

1.6288 

-240 

4.4969 

7.3202 

1.62783 

-180 

4.4948 

7.3201 

1.6286 

-260 

4.4968 

7.3198 

1.62778 

axis.  In  the  case  of  ice,  the  ratio  c/a  =  1.6293  (Megaw  1934).  In  case  of  balls,  the  corresponding  ratio  is 
c/a  =  2V2/V3  =  1.633. 

Owing  to  thermal  expansion,  parameters  a  and  c  (and  their  ratio)  change  with  temperature  (Table  3). 

When  “dumbbells”  are  packed  according  to  the  scheme ...  ABCABCABC  ...  the  structure  of  cubic  ice  1^ 
results,  which  is  the  closest  “relative”  of  ice  Ii,. 

The  syrmnetry  of  the  crystal  lattice,  consisting  of  oxygen  atoms,  is  characterized  the  space  group  Pbj/rrunc. 
In  this  notation,  “F*  indicates  that  the  space  lattice  is  primitive,  and  “63”  means  that  the  principal  symmetry 
axis  (c-axis)  is  a  hexagonal  screw  axis.  This  implies  that  if  an  array  of  oxygen  atoms  is  rotated  and  displaced 
simultaneously  along  that  axis  (in  a  screw-like  motirm),  such  that  the  displacement  is  c  per  revolution,  then  the 
first  coinciderice  will  occur  when  atoms  are  displaced  along  the  c-axis  by  c/2.  In  other  words,  the  symmeby  of 
the  c-axis  is  lower  than  six-fold.  In  the  case  of  a  pure  rotation  (without  translational  motion)  about  that  axis  by 
1/6  of  a  full  turn  (36076  =  60°),  the  coincidence  of  all  atcxns  will  not  take  place.  Otho-  symmetry  elements  in 
the  point  group  indicate  the  presence  of  planes  of  symmetry  (mirror,  m;  and  glide  plane,  c).  I  refo-  readers  who 
are  interested  in  a  detailed  description  of  these  elements  of  symmetry  to  one  of  the  textbooks,  Nye  (1989)  for 
instatKe.  Let  us  look  at  cases  when  these  elements  of  syrmnetry  are  somehow  related  to  the  symmetry  of 
particular  physical  prc^rties  of  ice. 

The  major  component  of  the  ice  crystal  lattice  is  a  water  molecule,  surrounded  by  four  neighboring  mole¬ 
cules  and  linked  to  them  by  hydrogen  bonds.  AtT-  0°C  the  distance  between  the  centers  of  the  closest  oxygen 
atoms  fgg  =2.76  A,  and  the  angles  between  hydrogen  bonds  are  very  close  to  the  angles  between  bisectors  of  a 
right  tetrahedron  -109°28'. 

A  unit  cell  of  I|,  structure  is  shown  in  Figure  13.  The  whole 
crystal  can  be  built  from  unit  cells  by  the  consecutive  process  of 
attaching  one  cell  to  another. 

Let  us  calculi  the  number  of  oxygen  atoms  per  unit  cell  of  I|, 
ice.  We  can  see  that  eight  atoms  are  located  in  the  vertices  of  the 
tetrahedron,  each  thus  belonging  to  eight  unit  cells  simultaneous¬ 
ly,  and  four  atoms  are  on  the  edges,  i.e.,  shared  among  four  unit 
cells.  Therefore,  there  are  8/8  +  4/4  +  2  =  4  oxygen  atoms  or  wa¬ 
ter  molecules  per  unit  cell. 

The  cocMdinate  system  defining  the  positions  of  all  oxygen  at¬ 
oms  as  well  as  crystallographic  planes  is  formed  by  lattice  vec¬ 
tors  coinciding  with  the  edges  of  the  unit  cell.  When  characteriz¬ 
ing  a  hexagonal  lattice,  it  is  more  convenient  to  use  four  rather 
than  three  unit  vectors,  as  is  illustrated  in  Figure  14.  Three  of  these 
vectors  lie  in  the  plane  perpendicular  to  the  hexagonal  c-axis,  and 

the  fourth  vector  coincides  with  the  c-axis.  Exan^les  of  coordi-  Figure  13.  Unit  cell  of  ice  /*. 
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nates  of  some  points  in  this  basis  are  shown  in  the  same  figure.  In 
the  MillCT-Bravais  system,  the  potential  redundancy  caused  by 
having  three  coplanar  lattice  vectois  is  resolved  by  ensuring  that 
the  first  three  coordinates  sum  to  zero. 

The  advantage  of  a  four-index  system  of  notation  is  that  all 
physically  equivalent  points  (in  the  case  of  hexagonal  symmetry) 
can  be  obtained  by  rearranging  the  first  three  indices.  For  instance, 
all  the  vertices  of  a  hexagon  have  coordinates  of  the  type  1/ 
3[1120],forexaii4>le,  1/3[2110],  1/3[1210]  and  soon. 

If  for  (hkil)  plane  either  1  =  0  or  h,k,i  all  =  0,  then  the  plane  is 
perpendicular  to  [hkU\  vector.  Thus,  ((X)01 )  denotes  the  plane  nor¬ 
mal  to  the  c-axis.  This  plane  is  called  the  basal  plane  and  is  the 
plane  of  the  closest  packing  in  I|,  structure.  Other  close-packed 
planes  are  prismatic  planes,  parallel  to  the  c-axis.  There  are  pris¬ 
matic  plaites  of  the  hrst  order  (1100)  and  of  the  second  order 
( 1 1 20).  Planes  that  form  an  angle  with  the  c-axis  and  ^90°) 

are  called  pyramidal. 

In  conclusion,  let  us  briefly  consider  the  similarities  atKl  dif¬ 
ferences  between  the  structures  of  ice  I|,  and  water.  Specifically, 
we  shall  try  to  understand  why  the  density  of  water  is  higher  than 
that  of  ice. 

Experimental  and  theoretical  studies  of  water  structure  car¬ 
ried  out  in  the  1970s  and  1980s  revealed  a  certain  random,  tetra¬ 
hedral  network  of  hydrogen  bonds  in  the  liquid  phase  (Malenkov 
1988).  The  average  distance  between  adjacent  water  molecules 
is  almost  the  same  as  in  the  case  of  ice,  but  the  average  number 
of  closest  neighbors  is  larger.  In  water  this  number  is  approxi¬ 
mately  4.4,  compared  to  4  in  ice  Ih.  This  dissimilarity  results  in  a 
larger  density  of  water  in  comparison  with  ice. 

There  is  a  widely  spread  opinion  that  the  lower  density  of  ice 
is  attributable  to  the  “openness”  of  its  structure,  in  which  there 
appear  to  exist  a  lot  of  “voids.”  Figure  15,  in  which  the  dimen¬ 
sions  of  atoms  are  given  in  the  same  scale  as  interatomic  distanc¬ 
es,  appears  to  support  this  contention. 


100011 


i  11010]  1(11201 

Figure  14.  Coordinate  sys¬ 
tem  for  hexagonal  lattices. 


Figure  15.  Same  view  of  the  ice  If,  lattice  as 
in  Figure  12b,  but  with  oxygen  atoms  rep¬ 
resented  by  balls  of  Van  der  Wools  radii. 


ARRANGEMENT  OF  PROTONS  AND 
CONFIGURATIONAL  ENTROPY  OF  ICE 

The  crystal  lattice  consisting  of  oxygen  atoms  considered  in  the  previous  section  does  not  exhaust  by  any 
means  the  description  of  the  ice  structure.  We  now  have  to  find  possible  positions  for  an  enormous  number  of 
protons  in  this  lattice.  This  question  is  not  trivial. 

Even  taking  into  account,  from  the  comparison  of  optical  spectra  of  water  vapor  and  ice,  that  water  molecules 
retain  their  properties  and  individuality  in  the  ice  lattice,  there  is  still  a  huge  number  of  possible  distributions  of 
protons  at  hydrogen  bonds  in  ice.  We  shall  try  to  exemplify  this.  Let  us  consider  an  ice  crystal  consisting  of  N 
water  molecules.  Since  the  number  of  hydrogen  bonds  in  the  crystal  is  Just  the  number  of  hydrogen  atoms,  then 
there  are  2N  hydrogen  bonds  in  this  crystal.  The  proton  on  each  bond  “must  know  its  master,”  i.e.,  it  must  be 
located  closer  to  either  one  oxygen  atom  or  the  other.* 

•  If  a  proton  were  positioned  in  the  middle  of  a  hydrogen  bond,  it  would  result  in  a  typical  ionic  crystal  with  a  much  higher 
melting  point  than  ice  (hundreds  of  degrees  Celsius),  huge  electroconductivity  of  the  liquid  phase  (water)  and  much  smaller 
dielectric  permittivity  t »  10,  whereas  in  the  case  of  ice  e  « 100. 
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Thus,  we  might  expect  (2)^  possible  arrangements  of  protons.  This,  of  course,  is  an  overestimated  value, 
since  there  will  be  many  incorrect  molecules  that  have  four,  three,  one  or  zero  protons  instead  of  two.  Still,  it  is 
clear  that  the  number  of  remaining  configurations  is  enormous. 

We  can  also  approach  the  problem  from  the  other  end,  focusing  our  attention  on  a  wato-  molecule.  If  we 
specify  the  positions  of  its  four  nearest  neighbors  in  space,  i.e.,  four  hydrogen  bonds,  there  are  still  six  possible 
orientations  of  the  molecule  left,  as  shown  in  Figure  16.  Since  the  total  number  of  molecules  is  N,  then  in  the  case 
of  their  independent  arrangement,  we  might  expect  (6)'^  various  combinations,  which  will  actually  differ  from 
each  other  only  by  proton  arrangement.  But  this  value  is  again  overestimated,  since  we  are  neglecting  the  order 
in  hydrogen  bonds,  and  among  them  there  would  be  wrong  bonds  with  two  protons  and  vdthout  protons  at  all. 

In  our  preceding  simple  discussion,  two  very  significant  principles  that  an  “ideal"  or  defectless  ice  structure 
must  obey  have  been  employed; 

1 .  There  are  two  protons  near  each  oxygen  atom. 

2.  There  is  a  single  proton  on  a  hydrogen  bond. 

The  fust  principle  states  the  electroneutrality  of  a  water  molecule  in  the  ice  crystal  lattice,  and  the  second  one 
takes  care  of  fulfilling  the  main  principle  of  a  hydrogen  bond,  the  presence  of  a  (single)  proton  at  it 

Both  principles  were  first  introduced  by  Bernal  and  Fowler  (1933),  whose  names  they  bear,  and  later  were 
more  seriously  grounded  and  considered  by  Pauling  (1935).  Sranetimes,  these  principles  are  just  referred  to  as 
the  “ice  rules.” 

L.  Pauling,  twice  a  Nobel  Prize  winner,  was  the  first  to  calculate  with  good  accuracy  the  numbor  of  possible 
configurations  of  protons  that  satisfies  the  ice  rules.  His  arguments  were  sinqile  and  elegant,  and  we  will  repro¬ 
duce  them. 

As  we  mentioned  above,  it  seems  at  first  glance  that  we  can  distribute  2N  protons  on  2N  hydrogen  bonds  in 
(2)^^  ways,  given  that  there  are  two  possible  positions  for  a  hydrogen  atom  on  each  bond. 

However,  if  we  consider  four  bonds  attached  to  an  oxygen  atom  (Fig.  16),  then  we  see  that,  instead  of  the 
expected  2^  =  16  configurations,  only  6  are  allowed.  This  means  that  for  each  lattice  node  a  “reduction  coeffi¬ 
cient”  of  6/16  must  be  used,  which  becomes  (6/I6)N  in  the  case  of  N  nodes.  Then,  the  total  number  of  configura¬ 
tions  W  will  be  given  by 

IV=(2)2Nx(6/16)N  =  (3/2)N  (3) 

The  value  of  W  is  still  not  absolutely  exact  for  two  reasons.  The  first  is  that  water  molecules  in  the  structure  of 
ice  Ih  are  incorporated  in  closed  rings,  the  smallest  of  which  comprises  six  molecules  (see  Fig.  12b,  for  instance). 
Therefore,  reorientation  of  a  single  molecule  in  the  ring  causes  reorientation  of  others,  if  the  ice  rules  are  not  to 
be  violated.  This  circumstance  is  extremely  uj  t ficult  to  take  into  account.  The  closest  approximation  was  report¬ 
ed  by  Nagle  (1966)  who  found 

W'=  (1.506851 0.000  15)N.  (4) 


This  deviates  by  less  than  0.5%  from  the  crude  result  (eq  3)  of  Pauling! 


Figure  16.  Six  possible  orientations  of  the  water 
molecule  in  the  ice  If,  lattice. 
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Another  circumstance  not  taken  into  account  by  Pauling  was  the  fact  that  real  ice  has  surfaces.  The  expiession 
for  IV  of  a  real  ice  crystal,  of  finite  dimensions  and  with  surfaces,  was  derived  by  Suzuki  ( 1966).  Since  the  ratio 
of  molecules  at  the  surface  to  molecules  in  the  bulk  tends  to  decrease  rapidly  with  increasing  volume,  in  most 
practical  situations  we  can  neglect  this  correction  to  W. 

The  above  analysis  implies  that  ail  allowed  configurations  have  the  same  energy  or  that  the  energy  difference 
is  small  in  comparison  with  k^T.  And  the  experimental  results  support  this  hypothesis.  But  why  are  the  Bemal- 
Fowler  configurations  almost  isoenergetic?  We  don't  know  tiiis  yet. 

We  have  considered  the  question  of  various  configurations  of  the  ice  protonic  subsystem  and  their  number 
with  this  degree  of  thoroughness  because  we  now  have  a  chance  to  com^dtt  the  predictioiis  of  this  hypothesis 
with  various  experimental  data. 

First,  if  ice  really  has  a  disoxlered  protonic  subsystem,  then  there  must  be  a  quite  large  entropy  associated 
with  this  disorder,  which  should  manifest  itself  in  the  thermal  properties  of  ke.  Such  entropy  we  will  refer  to  as 
configurational  entropy  S^.  Another  term  sometimes  used  for  this  quantity  is  zero-point  entropy. 

The  quantity  according  to  the  definition  of  entropy,  is 

=  (5) 

where  kg  is  Boltzman’s  constant. 

If  one  assumes  that  protonic  disorder  is  preserved  in  ice  at  the  lowest  temperatures  as  well,  then  can  be 
found  from  the  calorimetric  measurements  of  heat  Q  that  is  liberated  during  the  cooling  of  a  certain  amount  of 
dilute  water  vtqx)r  (diluteness  is  necessary  for  the  gas  to  be  considered  as  an  ideal  gas).  The  change  in  the  entropy 
of  the  system  is 


dS  =  SQ/T 


(6) 


and 


AS  = 


Ti 


5Q/T 


(7) 


Starting  an  experiment  with  water  vapor,  one  must  follow  its  consecutive  phases  of  condensation  into  the 
water  phase  and  crystallization,  measuring  the  heat  released  during  phase  transitions  and  that  released  because  of 
heat  capacity.  Therefore,  we  start  with  the  known  value  of  5  and  track  its  behavior  down  to  temperatures  close  to 
absolute  zero.  According  to  the  third  law  of  themradynamics 


limS=0  (8) 

7-*0 

i.e.,  one  might  expect  that  the  quantity  determined  in  the  experiment  described  above  goes  to  zero.  Such  mea¬ 
surements  have  been  performed  several  times  on  ice  (see,  e.g.,  Giauque  and  Stout  1936,  Flubacher  et  al.  1960, 
Sugisaki  et  al.  1968). 

Those  measurements,  instead  of  zero  as  expected  from  formula  8,  yielded 

HmS=  3.41  J  mol"' deg-'  (9) 

which  is  remarkably  close  to  the  value  predicted  from  fcxmulas  4  and  S 

Sc  =  kg  lnW= (3.41  ±  0.0008)  I  mol-'  deg-» .  (10) 

What  does  this  result  mean?  That  ice  violates  the  third  law  of  thermodynamics?  That  this  law  is  not  valid?  Of 
course  not.  This  confirms  that  a  disorder  in  the  protonic  subsystem,  and  the  configurational  entropy  associated 
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with  it,  exists!  As  we  shall  see  in  the  report  concenwd  with  other  forms  of  ice,  the  configurational  entropy  (or  a 
disrader  in  the  protonic  subsystem)  merely  does  not  have  enough  time  to  vanish  during  the  cooling  of  ice  for 
purely  kinetic  reasons.  The  reorientation  time  of  a  water  molecule  in  pure  ice  has  already  become  much  longa 
than  a  human’s  lifetime  at  liquid  nitrogen  temperature  (77  K).  However,  if  ice  is  slightly  doped  with  KOH, 
which  accelerates  the  process  of  molecular  reorientation  at  low  temperatures,  then  a  transition  occurs  at  a  tem¬ 
perature  of  about  72  K  to  a  partially  ordered  ice  XI,  which  liberates  of  most  of  the  entrc^y  (see,  e.g.,  Tajima  et  al. 
1982).  What  else  might  one  expect  to  discover  in  an  experiment,  if  such  a  disordered  or  “statistical”  structure  of 
the  (xotonic  subsystem  exists? 

According  to  the  model  described  above,  as  introduced  by  Pauling  and  which  is  usually  termed  a  statistical 
model,  a  proton  can  assume  one  of  two  positions  on  a  hydrogen  bond,  and  a  proton  switches  from  one  position  to 
the  other  many  times  per  second.  Since  the  orientation  of  the  electric  dipole  moment  of  a  water  molecule  is 
changed  because  of  this,  information  about  the  reorientation  time  of  water  molecules  may  be  provided  by  mea¬ 
surements  of  the  dielectric  permittivity  of  ice  as  a  function  of  the  frequency  of  the  external  electric  field.  If  the 
characteristic  time  that  it  takes  for  the  molecules  to  reorient  is  t,  then,  in  the  frequency  range  to  >=  lAi,  the  magni¬ 
tude  of  the  dielectric  permittivity  e  must  change  markedly.  This,  in  fact,  is  found  to  be  the  case  in  experiments, 
and  for  pure  ice  1|,  at  7= -Ky’C  the  value  of  t  =  5  x  10-5  seconds  was  obtained.  To  film  the  positions  of  protons, 
we  would  need  a  high  speed  camera,  operating  at  an  amazing  speed.  When  determining  proton  positions,  if 
experimental  techniques  are  used  that  require  more  time  than  T  for  the  acquisition  of  the  information  (electron  or 
neutron  diffiaction,  for  instance),  the  position  of  protons  is  averaged  and  at  each  point  (of  two)  on  a  hydrogen 
bond  we  find  a  “half-proton.” 

Such  measurements  conducted  by  several  research  teams  yielded  exactly  the  same  result.  It  was  Peterson  attd 
Levy  (I9S7)  who  were  the  first  to  apply  neutron  diffraction  to  determine  the  positions  of  hydrogen  nuclei  in 
monocrystals  of  ice  I|,.  Since,  as  a  result  of  the  interaction  of  the  magnetic  moment  of  a  neutron  with  the  random¬ 
ly  oriented  magnetic  momentum  of  a  proton,  a  rather  strong  scattering  occurs,  Peterson  and  Levy  conducted  their 

experiments  on  D2O  rather  than  on  ordinary  ice.  (The 
magnetic  moments  of  the  proton  and  the  neutron  com¬ 
pensate  each  other  to  a  large  degree  in  deuterium  nuclei.) 

An  arrangement  of  oxygen  and  hydrogen  nuclei  in  an 
ice  lattice  is  shown  in  Figure  17.  It  is  important  to  note 
that  O-H  and  H-H  distances  in  ordinary  and  heavy  ice 
can  be  slightly  different.  This  distinction  is  determined 
by  the  amplitudes  of  both  thermal  and  “zero-point”  os¬ 
cillations  differing  because  of  differences  in  the  mass  of 
H  and  D.  Owing  to  the  anharmonicity  of  the  oscillations, 
the  position  of  the  center  of  mass  depends  on  the  ampli¬ 
tude  and  appears  to  be  different  for  hydrogen  and  deute¬ 
rium. 

Extremely  accurate  measurements  of  proton  positions 
in  the  ice  I),  lattice  were  made  by  Kuhs  and  Lehmann 
(1981,  1983,  1987),  who  used  neutrons  with  a  small 
wavelength  X  =  0.7 107  A.  Distances  0-0,  O-H  and  H-H 
for  three  different  temperatures  are  summarized  in  Ta¬ 
ble  4.  Another  interesting  result  is  that  the  magnitudes 
of  proton  oscillations  were  determined  (see  Table  4). 
Owing  to  the  proton  having  the  smallest  mass  among  all 
atomic  nuclei,  the  magnitude  of  the  oscillations  turns  out 
to  be  very  significant,  even  at  very  low  temperatures. 
Therefore,  protons  can  be  viewed  as  being  “spread”  or 
“distributed”  over  a  region  with  dimensions  of  about  0.5 
A,  as  is  illustrated  in  Figure  18. 
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Figure  1 7.  "Half-proton  ”  picture  for  representing 
the  proton  arrangement  in  the  ice  I/,  structure  (after 
Peterson  and  Levy  1957). 
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Tibic  4  Moloadar  teasMtfy  ud  tfaennai  motfon  of  ice  The 
tinnili  ft  Mil  IT  ty  err  ilnt  tin  n  jHellinmiMi  T* iih  Thu 
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ki  ccflaiD  dinctloae  aliMi(  or  oillMitoaal  to  aik  o,  c  or  (ifler 

IfolMnllrlnn—  ITItl) 


60K 

123  K 

233  K 

Disttoccs  (A) 

o-ry 

2.750 

(I) 

2.751 

(I) 

2.759 

(2) 

o-or 

2.752 

(1) 

2.753 

(1) 

2.760 

(1) 

O-Hl 

1.004 

(2) 

1.008 

(2) 

'..008 

(4) 

0-H2 

1.005 

(1) 

1.004 

(1) 

1.004 

(2) 

Hl-Hl' 

0.741 

(3) 

0.735 

(4) 

0.744 

(7) 

m-HT 

0.741 

(1) 

0.745 

(2) 

0.753 

(3) 

Angles  (deg.) 


or-o-or 

109.33 

(2) 

109.34 

(2) 

109.35 

(4) 

or-o-or 

109.61 

(2) 

109.60 

(3) 

109.60 

(5) 

H1-0-H2 

109.21 

(11) 

10940 

(16) 

109.44 

(27) 

H2-0-H2' 

109.73 

(11) 

109.54 

(16) 

109.50 

(27) 

Rmsd  (A) 


0|c 

0.1179 

(5) 

0.1521 

(7) 

0.2069 

(14) 

Oic 

0.1183 

(9) 

0.1534 

(8) 

0.2096 

(10) 

Hl«a 

0.1337 

(16) 

0.1626 

(21) 

0.2115 

(41) 

HI  ic 

0.1784 

(6) 

0.1988 

(11) 

0.2409 

(28) 

HllOO' 

0.1358 

(10) 

0.1650 

(15) 

0.2109 

(31) 

H2|a 

0.1799 

(10) 

0.^ 

(14) 

0.2421 

(23) 

H2X<i 

0.1762 

(10) 

0.1984 

(15) 

0.2407 

(27) 

Employing  the  data  on  the  distribution  of  proton  density  along  a  hydrogen  bond,  Kuhs  and  Lehmann  estimat¬ 
ed  the  height  of  the  potential  barrier  separating  the  two  stable  positions  of  a  proton  on  a  bond.  They  obtained  for 
the  height  of  the  barriCT  the  very  small  value  of  0.012  eV,  which  can  be  easily  surmounted  by  a  proton,  even  at  T 
=  0  K,  due  to  quantum  tunneling.  However,  one  must  keep  in  mind  that  a  proton  cannot  make  such  a  transition 
from  one  well  to  another  alone,  because  the  first  ice  rule  would  be  violated.  Two  char^^  ions  would  appear, 
which  would  have  an  energy  of  about  1  eV,  much  larger  than  the  0.012-eV  barrio’  height.  Only  simultaneous 


Figure  18.  Probability  density  map  of  HI  in  the  (1010)  plme  as  described  by  har¬ 
monic  and  anharmonic  terms  up  to  fourth  order  (equidistant  contours,  zero  contour 
omitted).  The  small  negative  regions  are  due  to  minor  shortcomings  of  the  chosen 
mathematical  description.  With  the  limits  of  error,  these  probability  densities  are  fol¬ 
lowing  the  expectations  for  simple  temperature  dependence  assuming  Boltzmann  statis¬ 
tics  (after  Kuhs  and  Lehmann  1983). 
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(correlated)  motion  of  protons  along  closed  tnyectories  at  hydrogen  bonds  is 
allowed  by  the  ice  rules  without  the  formation  of  ions.  An  example  of  such  a 
motion  along  the  shortest  ring,  consisting  of  six  wata  molecules,  is  shown  in 
Figure  19.  The  question  of  the  possibility  of  collective  correlated  motion  of 
protons  in  ice  is  still  a  matter  of  controversy.  We  will  return  to  it  in  the  rqxxt 
concerned  with  point  defects  and  solitons. 

Shimaoka  (1960)  determined  proton  positions  in  ice  using  electron  difhac- 
tion.  Those  expoiments,  which  had  much  lower  accuracy  than  the  neutron 
diffraction  technique,  yielded  the  value  of  0.%  ±  0.03  A  for  the  O-H  distance. 

Finally,  Kuhs  and  Lehmann  (1987),  smdying  neutron  dif^action  in  ice  at 
very  low  (IS  K)  temperature,  showed  that,  because  of  the  “disorder”  in  the 
protonic  subsystem,  the  oxygen  lattice  is  also  disordered.  Namely,  the  exact 
positions  of  oxygen  atoms  deviate  randomly  (dqiending  tqxm  the  node)  from 
the  corresponding  positions  in  the  ideal  lattice  by  approximately  0.01  A. 


Figure  19.  Correlated  hops  of 
six  protons  along  the  hydro¬ 
gen  bonds. 


SPATIAL  NETWORK  OF  ELECTRONIC  DENSITY  OF 
ICE  AND  STRUCTURE  OF  ENERGY  BANDS 

One  can  say  that  ice  is  composed  of  the  atomic  nuclei  of  oxygen  atoms  and  hydrogen  atoms — protons  and 
electrons.  Now  we  have  to  consider  the  distribution  of  die  latter  in  the  ice  structure. 

From  a  methodological  point  of  view,  there  are  two  quite  different  apfMoaches  to  studying  the  electronic 
structure  of  ice.  The  first  one  is  based  on  the  use  of  certain  techniques  (X-ray  diffraction,  for  instance)  that  direc¬ 
tly  respond  to  the  magnitude  of  the  electron  density,  thereby  allowing  us  to  reconstruct  the  spatial  distribution  of 
electrons  in  the  ice  crystalline  structure. 

In  recent  years,  very  significant  results  have  been  achieved  in  this  area,  owing  to  both  advances  in  the  design 
of  X-ray  devices  and  improvement  of  techniques  for  analysis  of  data  (Goto  etal.  1990,  Sakata  etal.  1991, 1992). 
The  distribution  of  the  electronic  density  in  the  three  crystallographic  planes  of  hexagonal  ice  is  shown  in  Figure 
20,  adapted  from  the  paper  of  Sakata  et  al.  (1991),  who  used  for  the  accurate  structure  analysis  the  Maximum 
Entropy  Method  (MEM). 

Let  us  note,  in  the  first  place,  that  ice  is  permeated  by  a  rather  dense,  three-dimensional  network  of  electrons. 
This  network  “glues”  individual  water  molecules  into  a  rigid  crystalline  lattice,  w>'’ch  is  quite  distinct  from  the 
concept  of  ice  as  a  molecular  crystal  in  which  individual  molecules  interact  very  weakly. 

Information  about  the  distribution  of  electron  density,  such  as  that  represented  in  Figure  20,  is  extremely 
valuable  and  can  be  employed  to  derive  numerous  properties  and  parametos  of  the  material.  In  fact,  since  we 
already  know  the  arrangement  of  oxygen  and  hydrogen  atoms,  we  know  the  complete  distribution  of  electric 
charges  inside  the  ice  lattice.  This,  in  principle,  allows  us  to  calculate  the  energies  of  hydrogen  bonds,  the  struc¬ 
ture  of  filled  electronic  energy  bands,  elastic  constants,  etc.  However,  I  am  not  aware  of  conqiutations  of  tiiis 
kind,  which  are  based  upon  the  data  on  the  spatial  distribution  of  the  electron  density,  being  carried  out  for  ice. 

Nevertheless,  we  can  illustrate  the  relationship  between  the  distribution  of  electrons  in  space  and  their  energy 
characteristics  while  avoiding  the  use  of  sophisticated  theoretical  methods.  In  quantum  mechanics  there  is  a  rela¬ 
tionship  between  electron  density  and  the  minimal  kinetic  energy  of  electrons.  If  the  motion  of  an  electron  along 
the  x-axis  is  restricted  to  Ar,  then  the  uncertainty  in  its  momentum  according  to  Heisenbog’s  uncertainty 
principle,  cannot  be  less  than 

Apx>fr/Ax  (!0) 

where  h  is  Planck’s  constant.  The  corresponding  uncertainty  in  the  kinetic  energy  is 

AT  iAp^'ji^/2mg> hVlm^lAx)^  (11) 

where  mg  is  an  electron  mass. 
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a  MEM  map  o/(l  120) plane).  c.  MEM  map  of  (000,16)  plane. 


F igure  20.  Electron  density  distribution  of  ice  The  density  contours from  0  to  I A-^  are  shown  (ttfter  Sakata 

etal.  1991). 

In  the  case  of  motion  that  is  lestricted  in  three  dimensions  to  a  cube  with  side  Ax,  the  uncertainty  in  the  energy 
is  three  times  greater 

Ar= AT,  +  Afy  +  AT^  ^  3h2/2me(Ax)2.  (12) 

Formula  12  can  be  rewritten,  using  the  spatial  density  of  electrons  n  =  (Ax)-3 

AmH^rim^n^.  (13) 

The  width  of  an  individual  electronic  energy  band  in  a  solid  is  determined,  to  a  first  approximation,  in  terms  of 
the  difference  between  the  maximum  and  minimum  kinetic  energy.  Therefore,  the  result  of  formula  13  can  be 
interpreted  as  a  relationship  between  the  local  density  of  electrons  n  aitd  the  width  of  the  energy  band  (or  bands) 
to  which  the  electrons  belong. 

The  minimum  density  of  the  electrons  at  the  lines  of  hydrogen  bonds  (Fig.  20)  appears  to  be  n  =  3  x  1023 
cm-5,  fran  which 

ATS  3/|2/2me  n'i^  =  4.6  eV.  (14) 

Such  consideration  may  have  sense  only  if  the  electron  density  n  fills  the  whole  crystal,  so  that  we  can  use  free 
electron  gas  model  for  the  estimate. 

The  derived  value  (eq  14)  should  be  regarded  as  only  an  order  of  magnitude,  since  we  have  used  the  mass  of  a 
firee  electron  (in  a  vacuum)  m^  =  9.109  x  10-28  g,  and  genoaliy  one  should  use  the  effective  mass  of  the  electrons 
in  a  given  energy  band.  In  experimental  work  by  Shibagushi  et  al.  (1977),  the  total  width  of  the  first  three  filled 
bands  appears  to  be  equal  to  about  2  eV  and  the  maximum  difference  in  the  energies  for  the  electrons  in  those 
bands  is  about  6  eV.  Therefore,  the  crude  estimate  of  formula  14  turns  out  to  be  midway  between  these  two 
values.  Is  it  large  or  small?  It  is  well  known  that  in  isolated  stable  atoms  and  molecules  electrons  occupy  discrete 
energy  levels  of  zero  width.  As  M  molecules  are  brought  together,  each  of  these  levels  becomes  split  into 
closely  spaced  levels,  thereby  forming  an  energy  band.  When  bands  are  formed  witii  an  energy  width  of  several 
atomic  energies  (about  10  eV),  we  obviously  can  no  longer  speak  about  isolated  molecules  and  molecular  spec¬ 
tra  Consequently,  from  the  stanc^xnnt  of  the  electronic  structure,  ice  is  more  correctly  considered  a  rather  typical 
crystal  with  common  electrons  that  can  move  about  in  three  dimensions. 

The  second  methodological  approach  for  studies  of  the  electronic  structure  of  crystals  is  based  on  the  use  of 
experimental  methods  that  are  ciq)able  of  determining  tire  energy  characteristics  of  electrons:  the  position  of  the 
bottom  and  top  of  electronic  bands  with  respect  to  vacuum  level,  the  band  width,  the  position  of  tiie  Fermi  level 


Figure  21.  X-ray  photoelectron  spectra  of  ice.  Light  curve  shows  data 
of  water  vapor.  The  lower figure  shows  the  electronic  band  structure  of 
cubic  ice  calculated  by  Pastori-Parravicini  and  Resca  (1973)  (after 
Shibaguchi  et  al.  1977). 
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(the  chemical  potential  of  an  electron  in  a  ciystal),  the  width 
of  a  forbidden  band,  etc.  I  shall  list  such  experimental  ^h- 
niques  below,  along  with  examples  of  woric  in  which  those 
techniques  have  been  employed  for  ice. 

The  first  two  methods — X-ray  photoelectron  spectroscopy 
and  ultraviolet  photoelectron  spectroscopy — make  use  of  the 
external  photo  effect.  In  these  methods  a  beam  of  photons  with 
energies  greater  than  the  binding  energy  of  electrons  in  the 
material  under  investigation  is  absorbed  by  electrons  in  the 
crystal.  A  certain  fraction  of  these  electrons,  located  close  to 
the  surface,  is  knocked  out  into  a  vacuum  and  their  kinetic 
energy  is  measured.  Obviously,  since  the  energy  of  the  inci¬ 
dent  electrons  is  fixed,  the  kinetic  energy  of  the  electrons 
“knocked  out”  will  decrease  with  increasing  depth  of  energy. 
Thus,  the  kinetic  energy  spectrum  of  photoelectrons  contains 
infcmnation  about  how  deep  the  electron  energy  bands  are  and 
what  their  density  of  states  is. 

An  example  of  an  X-ray  photoelectron  spectrum  of  ice,  in 
which  at  least  four  peaks  are  clearly  seen  that  correspond  to 
the  four  most  shallow  bands  filled  with  electrons,  is  represented 
in  Figure  21 .  In  the  same  figure  the  spectrum  of  water  vapor 
and  the  results  of  theoretical  calculations  of  Pastori-Parravi- 
chini  and  Resca  (1973)  are  shown  for  comparison.  The  no¬ 
menclature  used  in  the  theoretical  diagram  will  be  explained 
later  in  this  section.  An  excellent  example  of  ultraviolet  pho¬ 
toelectron  spectra  is  represented  in  Figure  22. 

Another  efficient  source  of  information  about  the  structure 
of  the  electron  energy  bands  of  ice  is  the  spectrum  of  optical 
abstnption  in  the  range  of  vacuum  ultraviolet.  An  example  of 
such  a  spectrum  is  shown  in  Figure  23. 

Shibaguchi  et  al.  (1977)  employed  all  three  methods  char- 


Figure  22.  Photoelectron  energy  distribu¬ 
tion  curve  of  ice  for  (a)  hv  =  21.2  eV and  (b) 
hv  =  16.8  eV;  (c)  spectrum  of  ice  obtained 
from  curve  (a)  after  subtraction  of  the  back- 
ground  (after  Abbati  et  al.  1979). 


Figure  23.  Absorption  spectra  of  hexago¬ 
nal  ice  and  amorphous  ice  at  80  K  (after 
Sekietal.  1981). 
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Figure  24.  Energy  level  scheme  of  ice  revealed  by  X-ray  photo¬ 
electron  spectroscopy,  ultraviolet  photoelectron  spectroscopy 
and  vacuum  ultraviolet  absorption  spectra.  The  electronic  bands 
theoretically  calculated  by  Pastori-Parravicini  and  Resca  (1973) 
are  shown  as  a  rtference  (efter  Shibaguchi  et  oL  1977). 


acterized  above  to  establish  the  positions  of  energy  bands  of  ice.  The  results  of  their  analysis  are  shown  in  Figure 
24,  in  which  the  positions  of  energy  bands  are  given  relative  to  the  energy  of  a  free  electron  in  a  vacuum  (denoted 
as  vacuum  level).  According  to  this  schematic,  the  width  of  the  forbidden  band  of  ice  is  -  7.8  eV.  This  is  the 
distance  between  the  top  of  the  upper,  filled  band  and  the  bottom  of  the  lower,  empty  band 

Since  the  upper,  filled  band  has  no  vacant  states  and  the  width  of  the  forbidden  band  is  huge  compared  to  the 
thermal  energy  k^T  =  0.022  eV  (T  =  -10°C),  from  the  point  of  view  of  electronic  conductivity,  ice  must  be  an 
ideal  dielectric.  An  estimate  of  the  concentration  of  intrinsic  electrons  yields  n  S  lO-*^^  cm-^  at  r= -10^.  This 
means  that  in  the  whole  Antarctic  ice  cap  there  is  not  a  single  electron  that  could  conduct  current  and  that  would 
appear  to  be  due  to  a  thermoactivated  transition  from  the  valence  band  to  the  conductivity  band. 

At  present,  there  are  only  a  few  works  concerned  with  theoretical  calculations  of  the  energy  bands  of  ice.  The 
difficulties  of  such  calculations  arise  because  one  cannot  employ  the  powerful  methods  of  solid  state  physics 
based  on  the  use  of  symmetry  elements.  However,  the  structure  of  ice  lacks  even  a  translational  symmetry. 
Therefore,  the  first  calculations  of  the  electron  structure  of  ice  by  Pastori-Parravichini  and  Resca  (1971, 1973) 
were  performed  on  completely  polarized  ice  (which  has  a  maximum  electrical  dipole  moment).  In  addition, 
cubic  ice  possessing  higher  symmetry  than  hexagonal  ice  has  been  considered  in  their  calculations. 

Before  we  proceed  to  a  discussion  of  their  results,  let  us  remember  that  in  solid  state  physics  the  energy 
spectra  of  electrons  are  usually  pictured  as  the  curves  of  electron  energy  dependence  versus  their  momentum  p  or 
their  wave  number  k,  associated  with  the  momentum  via  the  relationship 

p  =  hk  (15) 

Those  dependencies  are  plotted  for  each  energy  band,  which  arises  frinn  discrete  levels  of  atoms  or  molecules 
due  to  their  interaction. 

Because  of  translational  symmetry,  the  dependence  of  energy  E(k)  appears  to  be  a  periodic  function,  with 
poiods  equal  to  vectors  of  the  reciprocal  lattice.  Because  of  that,  the  graphs  in  Figure  25  are  reduced  to  the  frrst 
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(a)  Energy  (eV) 


(b)  Energy  (eV) 


1?=  ^1.1.1)  T?=0  1?=  ^1.0,0)  1?=  ^0,1,0)  T?=0  T?=  ^1.0.0) 

Figure  25.  Valance  and  excited  bands  of  the  cubic-ice  model  structure  at  the 
points  shown  in  Figure  26  (after  Pastori  Parravicini  and  Resca  1973). 


Rrillouin  zone  (the  Wigner-Seitz  ceil  of  the  reciprocal  lattice). 
In  Figure  26,  the  Brillouin  zone  of  a  cubic  crystal  is  shown 
with  an  indication  of  cmain  points  at  its  boundary.  The  first 
few  energy  bands  of  a  cubic,  completely  polarized  ice,  com¬ 
puted  by  Pastori-Parravichini  and  Resca  (1973)  are  plotted  in 
Figure  2S.  The  model  used  by  them  has  been  refined  for  cal¬ 
culations  of  a  similar  structure  (Resca  and  Resta  1977).  These 
'‘y  authOTS  also  pointed  out  that  since  in  cubic  and  hexagonal  ice 
the  intermolecular  distaiKes  and  the  arrangement  of  four  near¬ 
est  neighbors  coincide,  the  major  essential  features  of  their  elec¬ 
tronic  energy  spectra  must  also  be  similar. 

MSignificant  progress  in  the  theoretical  calculations  of  the  elec¬ 
tronic  structure  of  ice  has  been  recently  achieved  by  Ryzhkin 
Figure  26.  Brillouin  zone,  symmetry  points  (*992),  who  developed  a  technique  that  takes  into  account  the 

and  symmetry  lines  of  the  cubic-ice  model  disordered  structure  of  hexagonal  ice.  It  turns  out  that  if 

structure  (after  Pastori-Parravicini  and  the  structure  obeys  Bemal-Fowler  rules,  then,  despite  the  ab- 
Resca  1973).  sence  of  complete  translational  symmetry,  the  energy  bands 

are  virtually  unchanged.  The  major  consequence  of  the  disor¬ 
der  is  that  it  will  give  rise  to  a  so-called  “localized”  band  in  the  middle  of  the  forbidden  band.  The  width  of  the 
forbidden  band  for  ice  I),,  according  to  Ryzhkin,  is  8.S  eV. 


LITERATURE  CITED 

Abbati,  I,,  L.  Brakovicb  and  B.  De  Mkhdis  (1979)  Investigating  ultraviolet  photoelectron  spectroscopy  of  ice. 
Solid  State  Communications,  29: 5 1 1-5 14. 

Bader,  R.F.W,  and  G  A.  Jones  ( 1  %3)  The  electron  density  disUibution  in  hydride  molecules.  I.  The  water  mole¬ 
cule.  Journal  of  Chemistry,  41:  586-606. 

Bernal,  J.D.  and  R.H.  Fowler  ( 1 933)  A  theory  of  water  and  ionic  solutions  with  particular  reference  to  hydrogen 


19 


and  hydroxyl  ion.  Journal  of  Chemical  Physics,  1:  515-548. 

Bertie,  J.£.,  L.D.  Calvert  and  E.  Whalley  (1%3)  Transformations  of  ice  II,  ice  EDI  and  ice  V  at  atmosfrfienc  pres¬ 
sure.  Journal  of  Chemical  Physics,  38:  840-846. 

Brill,  R.  and  A.  Tippe  (1967)  Gitterparameter  von  Eis  I  bie  Tiefen  Temperaturen.  Acta  Crystallographica,  23: 
343-345. 

Claydon,  C.R.,  GA.  Segal  and  HJS.  Taylor  (1971)  Theoretical  interpretation  of  the  optical  and  election  scat¬ 
tering  spectra  of  H2O.  Journal  of  Chemical  Physics,  54:  3799-3816. 

Fletcher,  N.H.  (1970)  The  Chemical  Physics  of  Ice.  Cambridge  University  Press. 

Flubacfaer,  P.,  AJ.  Leadbetter  and  JA.  Morrison  (1960)  Heat  capacity  of  ice  at  low  temperatures.  Journal  of 
Chemical  Physics.  33:  1751-1755. 

Giauque,  W.F.  and  J.W.  Stout  (1936)  The  entropy  of  water  and  the  third  law  of  thermodynamics.  The  heat 
capacity  of  ice  from  15  to  273°K.  Journal  of  the  American  Chemical  Society,  58:  1 144-1 150. 

Glaeser,  R.M.  and  C.A.  Coulson  (1965)  Multiple  moments  of  the  water  molecule.  Transactions  of  the  Faraday 
Society,  61:  389-391. 

Glen,  J.W.  (1974)  Physics  of  ice.  USA  Cold  Regions  Research  and  Engineering  Laboratory,  Monograph  n-C2a. 
Glen,  J.W.  (1975)  Mechanics  of  ice.  USA  Cold  Regions  Research  and  Engineering  Laboratory,  Monograph  II- 
C2b. 

Goto,  A.,  T.  Hondoh  and  S.  Mae  (1990)  The  electron  density  distribution  in  ice  Di  determined  by  single-crystal 
X-ray  diffnictometry.  Journal  of  Chemical  Physics,  93(2):  1412-1417. 

Hobbs,  P.V.  (1974)  Ice  Physics.  Oxford:  Clarendon  Press. 

Hoqjo,  G.,  N.  Kitamura,  K.  Shinuioka  and  K.  Mihama  (1956)  Low  temperature  specimen  method  for  electron 
diffraction  and  electron  microscopy.  Journal  of  the  Physical  Society  of  Japan,  11:  527-536. 

Kuhs,  W.F.  and  M.S.  Lehmaim  (1981)  Bond  length,  bond  angles  and  transition  barrier  in  ice  Di  by  neutron 
scattering.  Nature,  294: 432-434. 

Kuhs,  W.F.  and  M.S.  Lehmaim  ( 1983)  The  structure  of  ice  Di  by  neutron  diffraction.  Journal  of  Physical  Chem¬ 
istry.  87: 4312-4313. 

Kuhs,  W.F.  and  M.S.  Lehmann  ( 1987)  The  geometry  and  orientation  of  the  water  molecule  in  ice  Oi.  Journal  de 
Physique,  48(C1):  3-8. 

La  Placa,  SJ.  and  B.  Post  (1960)  Thermal  expansion  of  ice.  Acta  Crystallographica,  13:  503-505. 

Maeno,  N.  (1981)  Science  of  Ice.  Sapporo:  Hokkaido  University  Press  (in  Japanese). 

Malenkov,  G.G.  (1988)  Water.  Physical  Encyclopedia.  Moscow:  Soviet  Encyclopedia,  p.  294-97. 

Megaw,  H.D.  (1934)  Cell  dimensions  of  ordinary  and  “heavy”  ice.  Nature,  134(3397):  900-901 . 

Nagle,  J.F.  (1966)  Lattice  statistic  of  hydrogen-bonded  crystals  I.  The  residual  entropy  of  ice.  Journal  of  Mathe¬ 
matics  and  Physics,  J7:  1484-1491. 

Nye,  JJ.  (1989)  Physical  Properties  of  Crystals.  Oxford:  Oarendon  Press. 

Pastori-Parravicini,  G.  and  L.  Resca  (1971)  Symmetry  analysis  and  electronic  states  in  cubic  ice.  Journal  of 
P/iyjic5(C),4:L314-L317. 

Pastori-Parravadni,  G.  and  L.  Resca  (1973)  Electronic  states  and  optical  properties  in  cubic  ice.  Physical  Re¬ 
view  B,S:  3009-3023. 

Pauling,  L.  ( 1935)  The  structure  and  entropy  of  ice  and  other  crystals  with  some  randomness  of  atomic  arrange¬ 
ment.  Journal  of  the  American  Chemistry  Society,  57:  2680-2684. 

Peterson,  S.W.  and  H.A.  Levy  (1957)  A  single-crystal  neutron  diffraction  study  of  heavy  ice.  Acta  Crystallo¬ 
graphica,  10:  70-76. 

Resca,  L.  and  R.  Resta  (1977)  Energy  bands  in  cubic  ice.  Physica  Status  Solidi  (b),  81:  129-138. 

Ryzhidn,  LA.  (1992)  New  stages  in  understanding  ice  properties:  A  first  principle  approach.  In  Physics  and 
Chemistry  of  Ice  (Norikazu  Maeno  and  Takeo  Hondoh,  Ed.).  Hokkaido  University  Press,  p.  141-148. 

Sakata,  M.,  M.  Takata,  H.  Oshizumi,  A.  Goto  and  T.  Hondoh  (1991)  Electronic  structure  of  hydrogen  bond  in 
ice  Ih  obtained  by  the  maximum  entropy  method.  In  Summaries  of  International  Symposium  on  the  Physics  and 
Chemistry  of  Ice,  1-6  September  1991,  Sapporo,  Japan,  p.  1 1-12. 

Sakata,  M.,  M.  Takata,  H.  Oshizumi,  A.  Goto  and  T.  Hondoh  (1992)  Electronic  density  distribution  of  ice  Ih 


20 


obtained  by  the  maximum  entropy  method.  In  Physics  and  Chemistry  of  Ice  (Norikazu  Maeno  and  Takeo  Hon- 
doh,  Ed.).  Hokkaido  University  Press,  p.  62-68. 

Seki,  M.,  K.  Kobayashi  and  J.  Hakahara  (1981)  Optical  spectra  of  hexagonal  ice.  Journal  of  the  Physics  Society 
of  Japan,  50:  2643-2648. 

Shibaguchi,  T.,  H.  Onuld  and  R.  Onaka  (1977)  Electronic  structures  of  water  and  ice.  Journal  of  the  Physics 
Society  of  Japan,  42:  152-158. 

Shimaoka,  K.  (1960)  Electron  diffraction  study  of  ice.  Journal  of  the  Physics  Society  of  Japan,  15:  106-1 19. 
StiUinger,  F.N.  (1990)  Water.  Encyclopaedia  of  Physics  (R.G.  Lemer  and  L.G.  Trilf,  Ed.).  New  York:  VCH 
Publishers  Inc.,  p.  1353-1354. 

StiUinger,  F.N.  and  K.S.  Schweilzer  ( 1 983)  Ice  under  pressure:  Transition  to  symmetrical  hydrogen  bonds.  Jour¬ 
nal  of  Physical  Chemistry,  87: 4281-4287. 

Sugisaki,  M.,  H.  Suga  and  S.  Seki  (1968)  Calorimetric  study  of  the  glassy  state  IV.  Heat  capacities  of  glassy 
water  and  cubic  ice.  Bulletin  of  the  Chemical  Society  of  Japan,  41:  2591-2599. 

Suzuki,  Y.  (1966)  Disorder  entropy  of  ice.  Contributions  of  Institute  of  Low  Temperature  Science,  Hokkaido 
University,  A21:  1-44. 

Tsyinui,  J.,  T.  Matsuo  and  H.  Suga  ( 1982)  Phase  transition  in  KOH-doped  hexagonal  ice.  Nature,  299: 8 10-8 1 2. 
WhaUey,  E.  (1957)  The  difference  in  the  intermolecular  forces  of  H2O  and  D2O.  Transactions  of  the  Faraday 
Society,  S3:  1578-1585. 

WhaUey,  E.  (1976)  The  hydrogen  bond  in  Ice.  In  The  Hydrogen  Bond  (P.  Schuster,  G.  Zundel  and  C.  Sandorfy, 
Ed.).  North-Holland  Publishing  Company,  p.  1425-1470. 

WhaUey,  E.,  D.D.  Klug,  M.A.  Floriano,  E.C.  Svensson  and  V.F.  Sears  (1987)  Recent  work  on  high-density 
amorphous  ice.  Journal  de  Physique,  48(C1):  429-434. 

SELECTED  BIBLIOGRAPHY 

Chamberiain,  J.S.,  F.H.  Moore  and  N.H.  Fletcher  (1973)  Neutron-diffraction  study  of  H2O  ice  at  77  K.  In 
Physics  and  Chemistry  of  Ice  (E.  WhaUey,  S.J.  Jones  and  L.W.  Gold,  Ed.).  Ottawa:  Royal  Society  of  Canada,  p. 
283-284. 

Kobayashi,  K.  (1983)  Optical  spectra  and  electric  structure  of  ice.  Journal  of  Physical  Chemistry,  87:  4317- 
4321. 

Petrenko,  V.F.  (1993)  Electrical  properties  of  Ice.  USA  Cold  Regions  Research  and  Engineering  Laboratory, 
Special  Report  93-20. 


21 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


,  ataragci  nourparrMponM. 

inaMainkig  ttia  data  naadad.  and  oompMing  ant  lavMng  lha  oolacaan  of  Motmallatt  Sand  oonvnania  lagatdng  Mi  bufPan  aallmaia  or  any  (Mhaf  aapaa  ol  Ma  GolaGMan  ol  Mm 
lncludkioaiigoaallon«oriadiieinoMabUidan,totMaMng>anHaadqiiaflaraSai>itoaa,OaacloralatorMonnalionOpatallonaandRa|ioita  12lSJaaaiaonOa»<aHi|f«aay,8uaBi20«.Attaglan. 
VA  22202-4302,  and  Id  ttia  OfRoa  ol  Managaman  and  Budgal.  Papanmik  Raducdon  Piotael  (OTOi^iaa),  WaMnglan,  OC  20603. 


1.  AGENCY  USE  ONLY  (Lawre  biw*) 


2.  REPORT  DATE 

October  1991 


3.  REPORT  TYPE  AND  DATES  COVERED 


4.  TITLE  AND  SUBTITLE 

5.  FUNDiNG  NUMBERS 

Structure  of  Ordinary  Ice  I|, 

Part  I;  Ideal  Structure  of  Ice 

DAAL  03-91-0-0164 

6.  AUTHORS 

Victor  F.  Petrenko 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  A00RESS<ES) 

Thayer  School  of  Engineering  U.S.  Army  Cold  Regions  Research 

Dartmouth  College  and  Engineering  Laboratory 

Hanover,  New  Hampshire  03755  Hanover,  New  Hampshire  03755-1290 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

Special  Report  93-25 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Office 

Research  Triangle  Park 

Durham,  North  Carolina  27709 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

12a.  DtSTRIBUTION/AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is  unlimited. 

Available  from  NTIS,  Springfield,  Virginia  22161. 

13.  ABSTRACT  {Maximum  200  wonts) 

This  report  examines  the  defect-hnee  structure  of  ordinary  I|,.  reviewing  experimental  results  and  theoretical  models  of  the 
arrangement  of  oxygen  and  hydrogen  atoms  and  the  electronic  structure  of  ice.  Special  attention  is  given  to  the  results  of 
studies  from  the  last  two  decades.  Among  the  experimental  techniques  under  consideration  are  X-ray  diffraction,  optical 
absorption,  X-ray  photoelectron  spectroscopy  and  ultraviolet  optical  spectroscopy,  nuclear  magnetic  resonance  and  elec¬ 
tron  spin  resonance.  The  main  subjects  discussed  are  the  chemical  content  of  ice,  the  structure  and  inoperties  of  water 
molecules,  hydrogen  bonds,  phase  transitions,  protonic  disorder  and  electron  energy  spectrum  of  ice. 


14.  SUBJECT  TERMS 

Band  structure 
Hydrogen  bond 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


Ice 

Phase  diagram 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


Physics 

Structure 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


IS.  NUMBE^^  PAGES 


16.  PRICE  CODE 


20.  UMITATION  OF  ABSTRACT 


NSN  754001-280-5500 


Standard  Form  296  (Rav.  2-89) 

PtmoIWI  by  ANSI  8«1 239-18 
298-102 


